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Abstract 

We consider the production and two-photon decay of the CP-even Higgs bosons (h° and 
if ) of the Minimal Supersymmetric Standard Model (MSSM) at the Large Hadron Col- 
lider. We study in detail the dependence of the cross section on various parameters of the 
MSSM, especially the dependence on the mixing effects in the squark sector due to the 
Higgs bilinear parameter \x and the soft supersymmetry breaking parameter A. We find 
that the cross section for the production of these Higgs bosons has a significant dependence 
on the parameters which determine the chiral mixing in the squark sector. The cross sec- 
tion times the two-photon branching ratio of h° is of the order of 15-25 fb in much of the 
parameter space that remains after imposing the present experimental constraints. For the 
H° the two-photon branching ratio is only significant if the H° is light, but then the cross 
section times the branching ratio may exceed 200 fb. The QCD corrections due to quark 
loop contributions are known to increase the cross section by 50%. We find the dependence 
of the cross section on the gluon distribution function used to be rather insignificant. 
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1 Introduction 



It is well known that in the Minimal Supersymmetric Standard Model (MSSM) [[|], two 
Higgs doublets with opposite hypercharge are required in order to preserve supersymmetry. 
The physical Higgs boson spectrum in the MSSM consists of two CP-even neutral bosons 
h° and H°, a CP-odd neutral boson A and a pair of charged Higgs bosons H . The 
most important production mechanism for the neutral SUSY Higgs particles at the Large 
Hadron Collider (LHC) is the gluon fusion mechanism, pp — > gg — > h°, H°, A and the 
Higgs radiation off top and bottom quarks || . Except for the small range in the parameter 
space where the heavy neutral Higgs H° decays into a pair of Z bosons, the rare 77 decay 
mode, apart from rr decays, is a promising mode to detect the neutral Higgs particles, 
especially if b quark decays cannot be separated from the QCD background^]. 

This process was studied several years ago ||, and it was concluded that the lightest 
Higgs could be detected in this mode for sufficiently large values of the mass of the pseu- 
doscalar Higgs boson 3> mz- Similarly, the 77 channel is important for the discovery 
of H° for 50 GeV < < 150 GeV. Related studies have been presented in ref. |fl. 
While most of the emphasis in these earlier works was on the SSC, the discussion of the 
dependence of the cross section on various parameters is relevant. 

In this paper we present a fresh study of the hadronic production and subsequent 
two-photon decay of the CP-even Higgs bosons (h° and H°) of the MSSM, valid for the 
LHC energy of \/s = 14 TeV, and using gluon distribution functions based on recent HERA 
data |J, to reassess the feasibility of observing the CP-even Higgs bosons in this mode. 
As mentioned earlier, the gluon fusion mechanism is the dominant production mechanism 
of SUSY Higgs bosons in high-energy pp collisions throughout the entire Higgs mass range. 
We study the cross section for the production of the h° and H°, and their decays, taking 
into account all the parameters of the model. In particular, we take into account the mixing 
in the squark sector, the chiral mixing, which also affects the Higgs boson masses through 
appreciable radiative corrections. This was previously shown to lead to large corrections 

It should be kept in mind that if SUSY particles are light, decays to squarks and charginos, h° -► qq, 
fi° — > x + x~ , could be seen at LEP2. 
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to the rates 0. 

In the calculation of the production of the Higgs through gluon-gluon fusion, we 
include in the triangle graph loop all the squarks, as well as b and t quarks, the lightest 
quarks having a negligible coupling to the h°. On the other hand, in the calculation of 
decay of the Higgs to two photons, we include in addition to the above, all the sleptons, 
W ± , charginos and the charged Higgs boson. 

The Minimal Super symmetric Model contains several soft supersymmetry-breaking 
terms. We write the relevant soft terms in the Lagrangian as follows || 

[ V 2 m w cos p s/Zmwsmp J 

-M^Q - fhlvFu* - m 2 D d R] d R - M^tfjtfi - M^ 2 H ] 2 H 2 

+— |AA + A A j + — ^2 {A fc A fc + K k K k } . (1.1) 

Subscripts u (or U) and d (or D) refer generically to up and down-type quarks. The Higgs 
production cross section and the two-photon decay rate depend significantly on several of 
these parameters. 

Even without chiral mixing, two basic mass scales are required, those of squark and 
gaugino masses. The squark masses are determined by an SU(2)-doublet mass parameter, 
together with two SU(2)-singlet mass parameters, denoted in eq. (|1.1| ) by Mu, fhjj and 
ffiD, respectively. As is often done, we consider the case when all the SUSY-breaking 
squark masses are taken equal to a common mass parameter, which we denote by m 
(= Mu = mu = mo)- We shall consider the situation where this parameter is chosen to be 
150 GeV for the first two generations, and vary it over the values 150, 500 and 1000 GeV 
for the third generation. (Most of the plots presented will be for m = 500 GeV.) The 
physical squark masses will in general be different from these values, but they determine 
the relevant orders of magnitude. 

The contributions from the squark and Higgs sectors depend on the relative sign 
between the A and \i parameters, but not on their overall signs. I.e., we can choose 
A positive, but then we must study both negative and positive values of \x in order to 
cover all of parameter space. The Higgs sector depends on A and \x through the radiative 
corrections. 
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The contribution from the chargino sector depends on the relative sign between ji 
and Mi- The chargino contribution is independent of A. Thus, the h — > 77 decay rate is 
independent of the the over-all signs of A, n and M2, but depends on all the relative signs. 
To be more specific, we may choose M 2 positive, but then A and fi must both be allowed 
to take on negative and positive values in order to cover the full parameter space. In most 
of the parameter space though, the dependence on the chargino mass (and therefore also 
on the sign of M 2 ) is rather weak. In these regions it suffices to consider A positive. 

The signs of the off-diagonal terms in the squark mass matrices are given by the 
definition of A and /1, and also by the definition of the fermion masses. The sign chosen 
for the fermion mass terms show up in some of the couplings, in the fermion propagators 
and the spin sums. 

Having fi and A non-zero has implications for the Higgs masses as well. These will 
be presented in Sec. 2, together with constraints related to other masses, before we come 
to the study of the cross sections and decay rates in Sees. 3 and 4. 

2 Constraints on Parameter Space 

As discussed in the Introduction, there are various contributions to the production and 
decay of the lightest Higgs boson at the LHC, and, therefore, it is necessary to have a 
description of the full parameter space and the theoretical and experimental constraints 
on it before embarking on the calculation of the cross section and the decay. 

At the tree level, the CP-even neutral Higgs mass matrix is controlled by two 
parameters, which can be chosen to be tjia and tan/3. However, there are substantial 
radiative corrections [Q to the neutral Higgs masses which depend on, besides the top 
quark mass, the supersymmetric trilinear couplings (A u , A d ), the soft supersymmetry 
breaking masses (Mjj, rhu, m£>, etc.), the bilinear parameter \i in the superpotential, and 
tan/3 (= V2/V1, where t>2 and v\ are the vacuum expectation values of the two Higgs 
doublets of MSSM). More recent radiative corrections |]10| typically reduce the Higgs mass 
by 10-20 GeV. However, they are only valid when the squark masses are of the same order 
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of magnitude, and will therefore not be used in the present study. As long as the "loop 
particles" are far from threshold for real production, the cross section does not depend 
very strongly on the exact value of the Higgs mass. 

We shall assume that all the trilinear couplings are equal so that 



A u — A a 



A, 



(2.1) 



and take the top-quark mass to be 176 GeV []TTJ| in our numerical calculations. We vary the 
parameters which enter the neutral CP-even Higgs mass matrix in the following ranges: 



50 GeV < m A < 1000 GeV, 1.1 < tan/3 < 50.0, 
50 GeV < < 1000 GeV, < A < 1000 GeV. 



(2.2) 



Parts of the /z-tan f3 plane must be excluded because of the experimental constraints 
on the squark, chargino and h° masses. For low values of m, the lightest squark tends to be 
too light (below the most rigorous experimental bound, ~ 44.5 GeV |12|) or even unphysical 
(mass squared negative). The excluded region of the parameter space is indicated in fig. 1 
for m = 150 GeV, M 2 = 200 GeV, uia = 200 GeV and two values of the trilinear coupling 
.4. 

The allowed region decreases with increasing A, but the dependence on M 2 and 
is in this region rather weak. In order to have acceptable 6-squarks, /i and tan /3 must lie 
inside of the hyperbola-shaped curves. Similarly, in order to have acceptable t-squarks, 
the corners at large and small tan/3 must be excluded. 

The chargino masses are, at the tree level, given by the expression 
1 



m 2 ± 



2 
± 



(M 2 + /i 2 ) + m 2 w 

2 ■ 2 + m^y cos 2 2/3 + m 2 w (Mf + /i 2 



2/iM 2 sin2/3) 



1/2 



For fi = 0, the above expression simplifies to: 



m 2 ± 



\m\ +m 2 w ± 



^M 2 4 + m 2 w cos 2 2/3 + m 2 w M, 



nl/2 



(2.3) 



(2.4) 



For the case of \i = 0, we see that, for tan/3 3> 1, the lightest chargino becomes massless. 
Actually, small values of /i are unacceptable for all values of tan /3. The lower acceptable 



value for will depend on tan/3, but that dependence is rather weak. The region that 
is excluded due to the chargino being too light, increases with decreasing values of M 2 . 
We note that the radiative corrections to the chargino masses are small for most of the 
parameter space [I3|. We show in fig. 1 the contours in the /x-tan/3 plane outside of 



which the chargino has an acceptable mass (> 45 GeV) flifl . By the time the LHC starts 
operating, one would have searched for charginos with masses up to 90 GeV at LEP2. 
Contours relevant for LEP2 are also shown. 

For larger values of m, there is no problem with squark masses. However, then the 
radiative corrections to the Higgs masses get large, and correspondingly some regions of 
parameter space have to be excluded. This is illustrated in fig. 2, for in = 500 GeV. The 
corners at large values of |//| and tan (3 must be excluded since the h° mass there would fall 
below the experimental bound obtained at LEP flSfl . The extent of these forbidden corners 
grows rapidly as decreases below 0(150 GeV). They also increase with increasing values 
of A. 

The charged Higgs boson mass is given by 

m 2 H ± = + m\ + A (2.5) 
where A arises due to radiative corrections and is a complicated function of the parameters 



of the model [|16 



The radiative corrections to the charged Higgs mass are not, in general, as large as 
in the case of neutral Higgs bosons. This is due to an approximate global SU(2) x SU(2) 
symmetry [I7|, valid in the limit of no mixing. In certain regions of parameter space the 



radiative corrections can, however, be large. This is the case when the trilinear mixing 
parameter A is large, itla is small, and when furthermore tan /3 is large. We shall include 
the effects of non-zero A and /i in the calculation of the charged Higgs mass. The present 



experimental limit of order 40-45 GeV [I8| is not relevant, but presumably by the time 
the LHC starts operating, one will at LEP2 have searched for charged Higgs bosons with 
mass up to around 90 GeV. Even this bound does not appreciably restrict the parameter 
space as given in figs. 1 and 2. 

The neutralino mass matrix depends on four parameters. These are M 2 , M 1; fi and 
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tan (3. However, one can reduce the number of parameters by assuming that the MSSM is 
embedded in a grand unified theory so that the SUSY-breaking gaugino masses are equal 
to a common mass at the grand unified scale. At the electroweak scale, we then have |Tj| 

Mi = - tan 2 6 W M 2 (2.6) 

We shall assume this relation throughout in what follows. The neutralino masses enter the 
calculation through the total width of the Higgs boson. For the gaugino masses, we take 
M 2 to be 50, 200, or 1000 GeV. (Most plots will be for M 2 = 200 GeV.) The experimental 
constraint on the lightest neutralino mass rules out certain regions of the parameter space 
20| , but these depend on several parameters, and are therefore not reproduced in figs. 1 



and 2. They are generally correlated with the bounds on chargino masses |I3 



3 The lighter CP-even Higgs boson h° 

Let us first consider the cross section for 

pp^h°X (3.1) 

For M 2 = 200 GeV, rh = 500 GeV, and fi = 500 GeV, we show in fig. 3 this cross section 
for four values of A, the trilinear coupling parameter. The following features are rather 
striking: 

• The cross section decreases appreciably for large values of A. This is mainly due to 
an increase in the h° mass. 

• There are sharp edges at small values of tan/3, and also at small m^. The edges at 
small tan/3 are caused by the h° becoming light. At small values of rriA and large A, 
the couplings of h° to b quarks and r leptons become large, making the cross section 
very large in this region. 

For the same parameters as above, we show in fig. 4 the total decay rate, T(h° — > all) 
and the two-photon decay rate, T(h° — > 77). As opposed to fig. 3, here we only consider 
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two values of A, namely A = and A = 1000 GeV. The two-photon decay rate is seen 
to increase sharply at large values of A, but this does not result in a larger rate for the 
process 

pp -> h°X -> 77X (3.2) 

since the production cross section also decreases, as shown in fig. 3 (mostly due to an 
increase in the Higgs mass, m^o). In fig. 5 we show the cross section for the process ( |3.2| ). 
A characteristic feature of the cross section is that it is small at moderate values of uia, 
and then increases steadily with increasing tua, reaching asymptotically a plateau. This 
behaviour is caused by the contribution of the W to the triangle graph for h° — > 77. 
The h°WW coupling is proportional to sin(/3 — a), where a is defined in terms of masses 
(including radiative corrections) as 

„ / m 2 A — mi \ 7r , 

cos 2a = - cos 23 — =4 f- < a < 0. 3.3 

\m 2 H0 -m 2 h0 J 2 

For large tua, at fixed 8, all Higgs masses, except m^o, become large, so that h decouples. 
For large m^, we actually have sin(/3 — a) — > 1, which is why the cross section increases 
and reaches a plateau for large tua- 

In figures 6 and 7 we show contour plots of the cross section ( |3.2| ), for M 2 = 200 GeV, 
m = 500 GeV and \i equal to —500 and +500 GeV, respectively. For each case, four values 
of the trilinear chiral-mixing parameter A are considered, A = 0, 200 GeV, 500 GeV and 
1000 GeV (figure 7a is thus a different representation of figure 5). 

The /i-dependence of the cross section can for the case of M2 = 200 GeV and 
m = 500 GeV be described as follows. At moderate values (/i = ±200 GeV), there is 
not much difference between the cross section for positive and negative values of \i. The 
cross section has a significant dependence on ttia, being low at wia < O(300 GeV), then 
increasing steadily and reaching a plateau with increasing tua- The dependence on tan/? 
is rather weak. 

For increasing values of \fi\ (500 GeV, 1000 GeV) the change in the cross section 
is rather complex. This is basically caused by two phenomena: (1) At large values of 
the squarks become too light or unphysical, in analogy with the case (M 2 = 200 GeV, 
m = 150 GeV) shown in fig. 1. Hence, there are regions both at small and large values of 
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tan/3 where the cross section is not defined. (2) At large values of \fi\ and large values of 
tan/5 (all m^) the Higgs gets very light (due to radiative corrections). As a consequence, 
the cross section can get rather high, where not forbidden due to (1) above. 

The dependence of the cross section on M 2 and rh is described in table 1. For small 
values of m (~ 150 GeV), the possible ranges of tan j3, fi and A become severely restricted, 
in order to obtain physically acceptable squark masses. There is a significant increase in 
the cross section as m increases from 500 GeV to 1 TeV, to values of the order of 25-30 fb. 

For small values of M 2 (~ 50 GeV), the possible range in p, is restricted in order to 
obtain physically acceptable chargino masses. As M 2 increases beyond 200 GeV, there is 
little further change in the cross section. Details are given in table 1. 

The cross section has a modest dependence on the choice of gluon distribution func- 



tion used. For the plots shown here, we have used the recent GRV Set 3 distributions, 
which are the default of the PDFLIB. The BM Set 1 E2TI leads to an increase of the 



cross section by about 5-6%, whereas the MRS Set 29 (SO') |23| and CTEQ Set 24 (2pM) 



2~4| give reductions by 3-5% and 8-9%, respectively. These uncertainties are thus rather 
insignificant. 



4 The heavier CP-even Higgs boson H 

We shall here briefly consider the process 

pp^ H°X ^-rjX (4.1) 

which is of interest for small values of tua- 

The two-photon decay of the heavier CP-even Higgs proceeds dominantly through 
W loops, and is complementary to that of the lighter CP-even Higgs. It is only significant 
if rriA is small, hence m#o itself must also be light. At small itia the total H° decay rate 
is small and thus the branching ratio for it to go into two photons can be considerable. 
As a result, the cross section for the process ( |4.1| ) can at small values of uia reach values 
exceeding 200 fb. 
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Contour plots of the cross section are shown in figures 8 and 9, for four values of 
A, and for li = -500 GeV (fig. 8) and /i = 500 GeV (fig. 9). For li = -500 GeV there is 
a strong increase in the cross section with increasing values of A. For positive values of li, 
however, increasing values of A lead to a reduction of the cross section. At low values of 
rriA the Higgs mass ra#o is of the order of 110-140 GeV, and the h° mass is close to the 
experimental lower limit. 



5 Summary and concluding remarks 

We have discussed in some detail the cross section for producing the CP-even Higgs bosons 
at the LHC, in conjunction with their decay to two photons. Where the parameters lead 
to a physically acceptable phenomenology, the cross section multiplied by the two-photon 
branching ratio is for the lighter CP-even Higgs boson of the order of 20-30 fb. 

These numbers do not take into account QCD corrections. Such corrections have 
been evaluated for the quark-loop contribution, and lead to enhancements of the cross 
section of about 50% P5| . However, in the presence of chiral mixing the squark loops also 
contribute significantly. Since the QCD corrections for these are not available, we have 
decided it was more clean to simply leave out all higher-order QCD effects. One should of 
course keep in mind that they are very important. 

There is a modest increase of the cross section with increasing values of A (i.e., 
with increasing chiral mixing). This comes about as the result of two competing effects: 
with increasing A, the Higgs boson becomes more heavy, leading to lower production cross 
sections. This is however offset by a corresponding increase in the two-photon decay rate. 
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Table 1. Dependence of the cross section for pp — > h°X — > 77X on M 2 , m and 



M 2 


?7t — iou oev 

lightest squark too 
light for large values of 
tan/3 and increasing A 


m = 500 GeV 


in = 1000 GeV 


50 GeV 

lightest 
chargino too 
light for small 
values of tan (3 
and positive \i 


cross section significant 
only for small A, small 
tan /?, and small 
negative \i 


cross section lower than 
at M 2 = 200 GeV, 
especially at moderate, 
positive fi; for large \fi\ 
cross section significant 
only for narrow range 
in tan f3 


cross section larger 
than at rh = 500 GeV 


200 GeV 


cross section significant 
only for small A, small 
tan/3, and small 


"default" given in 
figs. 3-9 for 
H = ±500 GeV; 
complex dependence on 
tan p for larger \fj,\ 


cross section 
significantly larger than 
at m = 500 GeV, 
reaching well beyond 
25 fb for large range of 
fj,; less dependence on 
tan/3 than at 
m = 500 GeV, in 
particular for large A 


500 GeV 
1000 GeV 


very similar to 
M 2 = 200 GeV 


very similar to 
M 2 = 200 GeV 


very similar to 
M 2 = 200 GeV 
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Figure captions 

Fig. 1. Regions in the //-tan/3 plane which are ruled out by too light chargino (x ± ) 
and squark masses. The gaugino mass scale is M 2 = 200 GeV, m = 150 GeV, 
and rrtA = 200 GeV. At this low value of m, the squark masses are too light or 
unphysical in much of the /x-tan (3 plane. The hyperbola-like contours give regions 
that are excluded by the lightest b squark being below 45 GeV (solid) or 90 GeV 
(dashed). The more straight contours at large /i and small tan /3 similarly indicate 
regions that are excluded by the lightest t squark. In a) we consider the trilinear 
mixing parameter A = 0, whereas in b) we take A = 200 GeV. 

Fig. 2. Regions in the /i-tan/3 plane which are ruled out by too light chargino (x 1 * 1 ) and 
h° masses. Similar to fig. 1, but for m = 500 GeV. In a) we consider the trilinear 
mixing parameter A = 0, whereas in b) we take A = 1000 GeV. The solid 
(dashed) contours for small refer to the chargino mass m x ± = 45 (90) GeV. 
The unlabelled contours near the corners at large refer to regions where the 
h° mass would be below 45 GeV. 

Fig. 3. Cross section for pp — > h°X as a function of itia and tan/3 for M 2 = 200 GeV, 
m = 500 GeV, and fi = 500 GeV. Four values of A are considered: a) A = 0, 
b) A = 200 GeV, c) A = 500 GeV and d) A = 1000 GeV. 

Fig. 4. Total decay rate T(h° — > all) and two-photon decay rate T(h° — > 77), as functions 
of m A and tan f3 for M 2 = 200 GeV, m = 500 GeV, and fi = 500 GeV. Two values 
of A are considered: A = and A = 1000 GeV. 

Fig. 5. Cross section for pp — > h°X — > 77X as a function of itla and tan/3 for M 2 = 
200 GeV, m = 500 GeV, /1 = 500 GeV, and A = 0. 

Fig. 6. Dependence of the pp — > h° — > 77 cross section on and tan/3 for different 
values of the trilinear couplings A. Four values of A are considered: a) A = 0, 
b) A = 200 GeV, c) A = 500 GeV and d) A = 1000 GeV. Here M 2 = 200 GeV, 
m = 500 GeV, and /i = —500 GeV. The solid contours are at 15 fb, the long- 
dashed ones at 20 fb, and the short-dashed ones at 25 fb. 

Fig. 7. Dependence of the pp — > h° — > 77 cross section on and tan/3 for different 
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values of the trilinear couplings A. As fig. 6, except that \x = 500 GeV. 
Fig. 8. Dependence of the pp — * H° — * 77 cross section on m A and tan/3 for different 
values of the trilinear couplings A. Four values of A are considered: a) A = 0, 
b) A = 200 GeV, c) A = 500 GeV and d) A = 1000 GeV. Here M 2 = 200 GeV, 
m = 500 GeV, and fj, = -500 GeV. The contours are at 10 fb (solid), 20 fb, 
50 fb, 100 fb and 200 fb. 

Fig. 9. Dependence of the pp H° —> 77 cross section on m A and tan (3 for different 
values of the trilinear couplings A. As in fig. 8, except that fi = 500 GeV. 
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